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INTRODUCTION 


In  the  earliest  years  of  the  application  of  seismic  techniques  to 
the  Investigation  of  nuclear  explosions  and  their  discrimination  from 
earthquakes,  very  little  emphasis  was  placed  on  observations  made  at 
short  distances  from  the  seismic  source.  Recently,  however,  there  has 
arisen  the  possibility  of  negotiating  a  treaty  which  would  permit  the 
USA  to  install  seismic  stations  inside  the  USSR,  and  so  attention  is 
being  paid  to  the  analysis  of  measurements  made  at  regional  distances. 
With  the  techniques  which  have  been  developed  in  the  years  subsequent  to 
the  early  investigations  of  regional  seismology,  particularly  the 
advances  which  have  been  made  in  the  application  of  digital  computers  to 
a  variety  of  complex  problems,  it  may  now  be  possible  to  analyze 
quantitatively  many  of  the  complicated  features  of  regional  seismograms 
which  hampered  the  previous  studies. 

I 

If  regional  seismic  phases  such  as  P  ,  P  and  L  are  to  be  used  for 

n  f?  g 

event  location,  discrimination,  and/or  yield  estimation,  it  will  be 
necessary  to  understand  the  mechanism  responsible  for  their  generation 
and  the  effects  of  the  geological  medium  on  their  propagation.  For 
example,  the  rate  of  decay  of  the  signal  amplitude  with  distance  must  be 
known  in  order  to  determine  the  event  magnitude  from  the  observations, 
but  several  previous  studies  summarized  by  Pomeroy  et  al.  (1982)  have 
shown  that  the  decay  rates  of  regional  phases  are  strongly  dependent 
upon  the  particular  source-to-receiver  propagation  path  and  that  the 
amplitudes  are  also  strongly  dependent  upon  the  geological  structure  at 
the  observation  site  (Der  et  al.,  1981).  Path-  and  site-dependent 
effects  also  hinder  any  simple  interpretation  of  other  signal 
measurements,  such  as  the  coda  duration  or  the  phase  velocity,  in.  terms 
of  seismic  source  characteristics  such  as  depth.  If  regional  phases  are 
to  be  used  to  study  the  earthquakes  or  explosisons  which  generate  them, 
then,  it  will  be  necessary  to  evaluate  quantitatively  the  effects 
introduced  into  the  observed  records  by  the  propagation  paths  and  the 
receiver  sites.  The  effects  of  local  structure  at  the  receivers  have 
been  investigated  in  certain  previous  studies  (e.g.,  Barker  et  al., 
1980a;  Der  e_t  al_. ,  1980);  in  this  report  we  shall  examine  the  effects  of 
the  propagation  path.  On  account  of  its  structural  simplicity  and  its 
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relevance  to  the  geology  of  the  USSR,  we  shall  concentrate  on  an  earth 
model  of  a  tectonic  shield,  with  and  without  sedimentary  cover.  We 
shall  restrict  our  discussions  to  the  single  phase  L  ,  but  we  shall  also 
discuss  briefly  how  our  methodology  might  be  generalized  to  Incorporate 
other  regional  phases  as  well. 


BACKGROUND 


In  1952  Press  and  Ewing  identified  as  distinct  seismic  phases  two 
short-period  surface  waves  which  were  observed  to  propagate  across  North 
America,  The  phase  which  was  characterized  by  particle  motion  trans¬ 


verse  to  the  direction  of  propagation  was  designated  "L  ",  and  the  phase 
_  8 
which  exhibited  Rayleigh-type  particle  motion  was  designated  "R  ",  The 

8 

subscript  "g"  referred  to  the  granitic  layer  of  the  crust,  along  the 


base  of  which  the  two  phases  were  conjectured  to  propagate.  It  was  soon 
shown,  however,  that  these  surface  waves  are  observa t iona  1  ly  complex  and 
hence  that  the  simple  "granitic"  model  for  their  propagation  is 


inadequate.  In  particular,  Bath  (1954)  divided  the  L  wavetrain  into 

8 

two  distinct  arrivals,  which  he  denoted  as  "L  1"  and  "L  2",  and  Bath 

8  8 

(1954)  and  Gutenberg  (1955)  showed  that  the  L  particle  motion  is  not 

8 

simply  transverse  but  that  it  has  large  vertical  and  radial  components 


as  well.  These  investigations  proposed  that  the  observations  should  be 


explained  by  considering  L  and  R  to  be  channel  waves  in  the  crust  or 

8  8 


upper  mantle,  i.e.,  waves  which  propagate  along  a  low-velocity  zone 
(LVZ)  waveguide  capped  by  a  higher-velocity  lid,  A  major  step  towards 


understanding  the  true  nature  of  L  was  taken  by  Oli  er  and  Ewing  (1957; 

8 

1958a, b),  who  noted  that  the  group  velocities  of  prominent  maxima  in  the 


L  wavetrain  correspond  to  stationary  points  in  the  group  velocity- 

o 


versus-f requency  dispersion  curves  for  higher-mode  Love  waves.  The 


distinct  arrivals  such  as  L  1  and  L  2  were  therefore  identified  as  Airy 

8  8  7 


phases,  and  their  amplitudes  and  rates  of  decay  could  be  predicted  using 
well-known  formulas  derived  by  means  of  the  stationary  phase  approxi- 

o 

mation  (e.g,,  Bath,  1968,  Chapter  3). 


The  identification  of  the  transverse  component  of  L^  as  the  super¬ 
position  of  higher-mode  Love  waves  has  permitted  the  observed  group 
velocities  of  that  regional  phase  to  be  used  as  a  tool  in  the  computa¬ 
tion  of  earth  models,  since  the  values  of  the  group  velocity  minima  in 
the  dispersion  curves  are  sensitive  to  the  details  of  the  layered 
structure  which  is  responsible  for  the  Love-wave  propagation.  Kovach 
and  Anderson  (1964)  calculated  the  group  velocities  for  low-frequency  L 
waves  in  order  to  test  various  assumed  values  for  the  thickness  and  the 
shear-wave  velocity  of  the  LVZ  in  the  mantle.  For  the  low-frequency 


provide  a  plausible  explanation  of  the  failure  of  L  to  propagate;  across 

long  oceanic  paths.  Although  the  shield  model  which  we  shall  consider 

in  this  report  contains  no  crustal  LVZ  and  although  L  is  unaffected  by 

8 

the  mantle  LVZ  at  the  high  frequencies  which  we  shall  consider,  we  shall 
find  this  same  method  of  analysis  useful  for  explaining  the  differences 
in  L  propagation  which  are  introduced  when  an  attenuating  sedimentary 

o 

layer  is  added  to  the  shield  structure. 


An  important  development  in  the  analysis  of  the  phase  L  was  the 

8 

generation  of  synthetic  seismograms  by  Knopoff  et  al.  (1973).  These 
investigators  used  the  phase  ve 1 oc i ty- vers  us- f requency  relations  and  the 
stress  and  particle  d  isp  laceme n t-versus-d e p th  relations  described 
previously,  and  they  combined  them  with  the  excitation  coefficients  for 
the  various  higher  modes  wh ich  were  compu  ted  by  assuming  various 
double-couple  source  mechanisms,  The  Love-wave  higher  modes  were  Bummed 
and  transformed  from  the  frequency  domain  into  the  time  domain,  result¬ 
ing  in  synthetic  seismograms  which  exhibited  many  of  the  observed  fea¬ 
tures  of  L  .  The  synthetics  were  a  convincing  demonstration  of  the 
8 

identification  of  L  as  the  sum  of  higher-mode  Love  waves,  an  identi- 

o 

fication  which  had  previously  been  surmised  simply  on  the  basis  of  the 
group-velocity  minima,  since  the  synthetics  displayed  not  only  the 
proper  group  velocities  but  also  the  proper  amplitudes  and  frequency 
content.  Knopoff  et  al.  (1973)  also  showed  by  means  of  these  synthetic 

O 

seismograms  that,  contrary  to  the  original  conjectures  of  Bath  (1954) 
and  Gutenberg  (1955)  but  in  accordance  with  the  results  of  Kovach  and 
Anderson  (1964),  no  crustal  LVZ  is  necessary  for  the  propagation  of  L  . 
Realistic  L  seismograms  can  be  synthesized  by  assuming  a  crustal  model 
composed  of  layers  for  which  the  shear-wave  velocity  increases 
monotonical ly  with  depth.  In  this  report  we  shall  verify  that  this 
result  holds  at  higher  frequencies  than  the  cutoff  (f  ■  1.0  Hz)  which 
was  apparently  used  by  Knopoff  et  al.  (1973). 

Computation  of  higher-mode  Love-wave  group-velocity  minima  and  L 

8 

synthetic  seismograms  was  also  used  by  Knopoff  et  al,  (1974)  to  deter¬ 
mine  shear-wave  velocity  profiles  for  different  types  of  continental 
crustal  structures.  Because  the  synthetics  were  generated  for  the 
long-period  WWSSN  instrument,  lower  frequencies  and  deeper  earth  models 
were  used  in  that  study  than  will  be  considered  in  this  report,  but  many 
of  their  results  are  important  for  understanding  the  characteristics  of 
L  which  are  observed  on  short-period  instruments  as  well.  In  parti- 

O 

cular,  Knopoff  et  al,  (1974)  found  that  both  the  minimum  group  velocity 
and  the  period  of  the  lowet t- f requency  Airy  phase  for  each  Love-wave 
mode  increase  monotonica 1 ly  with  crustal  thickness.  By  scaling  the 
crus  tal  thickness  from  a  continental  to  an  oceanic  earth  model, 
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Knopoff  et  al.  (1979)  showed  why  does  not  propagate  through  even 
short  paths  of  oceanic  crust  (ss  opposed  to  the  upper  mantle,  which  wss 
investigated  by  Schwab  and  Knopoff  (1971)).  The  absence  of  oceanic  L 

g 

can  be  attributed  Lo  the  scaling  with  crusta!  thickness  of  the  dis¬ 
persion  curves,  surface-wave  scattering  due  to  lateral  variations  in 
oceanic  sediment  thickness,  and  anelastic  attenuation  in  the  sediments. 
These  same  features  will  also  be  important,  although  to  s  lesser  extent, 
in  our  examination  of  L^  propagation  through  attenuating  sedi¬ 
ments  in  a  continental  earth  model. 

All  of  the  analyses  of  L  which  have  been  mentioned  thus  far  have 

g 

considered  particle  motion  in  only  the  trsnsverse  direction,  This 
restriction  was  brought  about  by  the  modeling  of  L  as  the  superposition 
of  msny  Love-wave  modes,  all  of  which  have,  of  course,  only  a  transverse 

O 

component  of  motion.  It  was  shown  by  Bath  (1954)  and  Gutenberg  (1955), 

however,  that  L  has  a  significant  vertical  component  of  motion.  A 

quantitative  explanation  of  both  this  and  vertical  and  the  radial 

components  of  L^  was  provided  by  Panza  and  Calcagnile  (1975),  who 

demonstrated  that  these  components  are  the  superposition  of  higher-mode 

Rayleigh  waves.  They  also  demonstrated  that  the  phase  R  is  the 

8 

fundamental  Rayleigh  mode  and  that  no  LVZ  in  either  the  crust  or  the 

upper  mantle  is  necessary  for  the  propagation  of  either  R  or  the 

g 

Rayleigh-type  of  L^  motion,  Although  we  shall  confine  our  synthetic 
seismograms  to  the  superposition  of  Love-wave  modes  and  hence  to  the 
transverse  component  only,  we  shall  show  the  necessity  of  including  the 
Rayleigh  modes  in  order  to  reproduce  observed  signals.  If  the  phase  L 

8 

is  to  be  used  as  a  seismic  discriminant,  it  will  be  important  to 
understand  the  partitioning  of  surface-wave  energy  between  SH-type  and 
P-SV-type  motion. 


In  this  report  we  have  undertaken  to  utilize  the  modal  superposition 
technique  in  order  to  answer  some  questions  about  L^  propagation  which 
ought  to  be  understood  if  L^  is  to  be  used  in  a  program  of  seismic 
discrimination  and  yield  determination  at  regional  distances.  In  parti¬ 
cular,  it  is  well  known  that  L^  not  only  fails  to  propagate  through  the 
ocean  but  also  attenuates  rapidly  in  certain  continental  regions.  It  is 

a  1  so  known  that  the  amp litude  of  L  ,  as  well  as  its  group  velocity,  is 

g 


strongly  influenced  by  surficial  geology.  In  order  to  interpret 

observations  of  7,  in  terms  of  the  source  mechanism,  as  well  as  to 

8 

choose  the  best  sites  for  locating  future  seismic  stations,  it  is 
important  to  have  a  quantitative  estimate  of  these  path  effects.  We 
should  also  like  to  know,  for  a  given  source-to-receiver  earth  model, 
how  large  a  variability  in  the  characteristics  of  L  can  be  anticipated 

s 

from  earthquakes  of  different  source  mechanisms  and  of  different  depth. 
It  is  especially  important  to  know  whether  certain  earthquakes  have  some 
L  signal  characteristics  which  can  be  used  as  diagnostics  in  discrim- 

o 

ination  studies,  Our  study  of  L  synthetic  seismograms  has  been  aimed 

8 

at  investigating  these  questions. 
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SKI. SMOG RAM  SYNTHESIS 


The  theory  of  surface-wave  propagation  through  a  stratified  medium 
i®  addressed  in  a  host  of  previous  studies;  herein  we  shall  discuss  only 
a  few  points  of  the  theory  which  are  relevant  to  the  generation  of  the 
synthetic  seismograms.  Our  synthetics  were  calculated  using  formulas 
for  far-field  surface  waves  given  by  Mendiguren  (1977),  which  we  shall 
not  reproduce  herein.  The  actual  computation  of  the  quantities  which 
appear  in  these  formulas,  however,  was  performed  using  a  computer 
program  based  on  the  method  of  Harkrider  (1970).  We  are  able  to  make 
this  substitution  since  the  results  obtained  by  both  methods  are 
eqivalent,  although  the  method  of  Mendiguren  (1977)  permits  the  use  of  a 
more  general  seismic  source  function  than  does  the  method  of  Harkrider 
(1970),  which  considers  only  certain  specific  source  functions.  Since 
all  of  our  synthetic  seismograms  will  be  based  upon  one  of  these  special 
cases,  namely  a  double-couple  source  mechanism,  the  two  formulations  may 
be  considered  identical  except  for  differences  in  notation  and  sign 
conventions . 

In  order  to  generate  synthetic  surface-wave  seismograms,  the 
following  quantities  must  be  computed: 

1)  the  seismic  source  function.  Since  our  synthetics  are  meant  to 
be  valid  for  the  far-field  displacement  only,  we  assume  that  the  seismic 
source  Is  of  inf inltessimal  spatial  extent  and  thus  may  be  represented 
by  a  second-order  moment  tensor,  .  For  the  special  case  of  a 
double-couple  mechanism,  the  nine  terms  of  this  moment  tensor  (only  six 
of  which  are  distinct)  can  be  expressed  In  terras  of  the  strike  and  dip 
of  the  fault  plane  and  the  rake  of  the  slip  vector.  We  should  point  out 
that  if  an  isotropic  explosion  rather  than  a  double  couple  were  used  as 
the  source  mechanism,  the  moment  tensor  would  become  simply 
Mi1  "  Mo  ^ij»  in  whlch  case  there  is  no  shear  displacement  at  the 
source  and  all  Love-wave  motion  would  therefore  disappear.  We  are  thus 
unable  to  model  the  transverse  component  of  L  from  an  explosion,  since 
It  is  the  departures  from  the  idealized  case  which  give  rise  to  the 
transverse  component  of  motion,  and  we  shall  restrict  our  analyses  to 
idealized  earthquakes  instead.  We  shall  assume  that  the  seismic  source 


Is  a  temporal  step  function,  although  we  can  generate  the  synthetics 
equally  well  using  some  other  source  time  function  if  we  wish  to  do  so. 

2 )  the  phase  velocity  of  the  surface  waves  as  a  function  of 
frequency .  For  a  given  earth  model  consisting  of  a  vertical  sequence  of 
constant-thickness,  flat  strata  in  which  the  density  and  the 
compresslonal-  and  shear-wave  velocities  are  specified,  the  computation 
of  the  Love-wave  and  Rayleigh-wave  phase  velocities  for  any  frequency  is 
an  eigenvalue  program  (the  "period  equation")  which  may  be  solved  by  the 
well-known  matrix  method  expounded  by  Haskell  (1953)  and  formalized  by 
Gilbert  and  Backus  (1966).  We  will  find  a  sequence  of  roots  as  the 
solution  to  this  eigenvalue  problem,  corresponding  to  the  sequence  of 
fundamental  and  higher  modes. 

3)  the  transfer  function  of  the  earth  model.  This  term,  denoted 
for  the  Love  and  Rayleigh  waves  as  A^  and  A^  by  Harkrlder  (1970)  and  as 
1/2C  0  I  ^  and  1/2C^U^T^  by  Mendiguren  (1977),  is  inversely 
proportional  to  the  integral  over  depth  of  the  kinetic  energy  density  of 
the  surface  wave.  This  term  shows  the  excitation  of  motion  at  the 
surface  by  a  wave  with  a  given  frequency  and  a  given  mode  number. 

4 )  the  stress  and  particle  displacement  as  a  function  of  both  fre¬ 
quency  and  depth.  For  the  Love  wave  this  involves  only  the  transverse 
direction,  but  for  the  Rayleigh  wave  both  the  vertical  and  radial 
components  must  be  calculated.  The  stress  and  the  displacement  may  be 
found  by  substituting  back  into  the  period  equation  each  root  which  was 
found  for  that  equation,  resulting  in  an  eigenfunction  for  every 
frequency  and  each  mode.  This  computation  involves  a  product  of 
matrices,  one  matrix  for  each  layer  in  the  earth  model,  which 
"propagate"  the  surface  boundary  conditions  downward  to  the  focal  depth 
of  the  source.  This  product  of  "propagator"  matrices  may  be  numerically 
unwieldy,  since  round-off  errors  will  be  compound'3-  by  each  successive 
multiplication  and  thus  will  increase  with  depth.  We  have  therefore 
found  it  necessary  to  impose  a_  priori  upon  the  computation  the 
constraint  that  the  eigenfunctions  in  the  bottom  layer  decay 
exponentially.  The  error  in  the  multiplication  of  the  propagator 
matrices  will  of  course  be  increased  greatly  by  inaccuracy  in  the 
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calculation  of  the  roots  of  the  period  equation.  The  calculation  of 
these  roots  is  in  fact  an  error-prone  process,  since  the  eigenvalue 
equation  is  based  on  the  vanishing  of  a  determinant  which  is  subject  to 
significant  numerical  imprecision  due  to  round-off.  For  the  Love  waves, 
the  propagator  matrices  describe  the  stress  and  displacement  in  only  one 
component  and  thus  are  2  x  2  matrices;  for  the  Rayleigh  waves,  there  are 
tve  components  and  so  the  matrices  are  4x4.  causing  the  determinantal 


equation  to  be  more  difficult  to  solve  accurately  than  is  the  case  for 

the  Love  waves,  particularly  at  high  frequencies.  The  imprecision  in 

the  eigenvalues  leads  to  eigenfunctions  which  are  sometimes  ill  behaved 

with  depth.  It  is  therefore  easier  to  synthesize  the  Love-type  motion 

of  L  than  the  Rayleigh-type  motion. 

8 


5)  the  azimuth  and  distance  from  the  seismic  source  to  the  point 
of  observation.  The  azimuth,  measured  with  respect  to  the  same 
reference  frame  as  are  the  source  function  coefficients  ,  gives  the 
location  of  the  observer  on  the  radiation  pattern  of  the  source 
mechanism  for  Love  and  for  Rayleigh  waves.  We  shall  make  the  distance 
dependence  of  the  amplitudes  of  our  synthetic  seismograms  slightly  more 
realistic  by  replacing  Mendiguren's  (1977)  assumption  of  cylindrical 
wavefront  spreading  with  a  term  appropriate  to  the  spreading  of  a 
wavefront  on  the  surface  of  a  sphere: 


amplitude 


o  (R  sin 


(rrt)) 


-1/2 


where  r  Is  the  source-to-receiver  distance  and  R  is  the  radius  of  the 

o 

earth.  We  reserve  for  later  a  discussion  of  how  the  decay  of  amplitude 
with  distance  Is  affected  by  anelastic  attenuation. 

On  account  of  its  prevalence  in  many  different  regions  of  the 
world,  Including  the  USSR,  the  earth  model  which  we  shall  use  for 
studying  L  propagation  is  that  of  a  tectonic  shield.  We  shall  take  our 

O 

model  from  the  study  of  the  Canadian  shield  by  Brune  and  Dorman  (1963). 
The  use  of  this  particular  shield  model  may  facilitate  the  comparison  of 
our  synthetics  with  observed  seismograms  recorded  at  North  American 
stations,  and  we  anticipate  that  the  synthetics  generated  with  this 
moael  will  have  many  salient  features  in  common  with  those  which  would 
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bo  generated  using  shield  models  appropriate  to  other  continents.  The 

model  of  Brune  and  Dorman  (1963)  is  reproduced  in  Table  I.  We  note  that 

since  our  Lg  synthetics  are  confined  to  Love-wave  motion,  we  shall  not 

use  the  tabulated  values  of  the  P-wave  velocity.  The  shear-wave 

velocity,  density,  and  modulus  of  rigidity  are  illustrated  as  functions 

of  depth  in  Figures  1-3.  It  is  shown  in  those  figures  that,  as  a 

computational  aid,  we  have  subdivided  the  bottom  layer  of  the  crust  so 

that  the  transition  to  the  mantle  is  slightly  less  abrupt  than  is  the 

case  in  Table  I.  This  subdivision  of  the  transition  zone  will  improve 

the  accuracy  of  the  computations  at  high  frequencies,  since  Imprecision 

can  result  if  the  exponentially  attenuating  layer  in  which  the  particle 

motion  dies  out  is  allowed  to  be  many  wavelengths  thick.  By  making  the 

base  of  the  crust  into  several  layers,  this  problem  is  avoided,  since 

the  program  Ignores  the  sublayers  below  the  depth  at  which  (for  a  given 

frequency)  the  particle  motion  vanishes.  We  see  that  there  is  no 

crustal  LVZ  for  this  shield  model.  The  mantle  LVZ,  which  lies  at  a 

depth  of  greater  than  100  km,  will  have  little  effect  upon  L  at  the 

frequencies  which  we  shall  consider  in  this  report.  For  this  hard-rock 

model,  a  value  of  1000  has  been  chosen  for  shear  wave  figure  of  merit  Q, 

independent  of  both  frequency  and  depth.  Such  an  assumption  represents 

the  limiting  case  of  only  slight  attenuation  of  l  . 

g 


Figure  4  shows  how  the  phase  velocity  of  each  the  first  twenty-six 
Love-wave  modes  behaves  as  a  function  of  frequency  from  0.02  to  2.50  Hz, 
For  the  sake  of  visual  clarity,  the  odd-numbered  modes  (counting  the 
fundamental  as  mode  1)  are  shown  as  solid  lines  and  the  even-numbered 
modes  as  shown  as  dashed  lines.  Phase  velocities  greater  than  4.51 
km/sec  are  not  considered,  because  this  cut-off  corresponds  to  the 
shear-wave  velocity  in  the  mantle  LVZ  (cf.  Table  I).  The  Love-wave 
modes  osculate  at  this  phase  velocity,  so  the  separate  curves  shown  in 
Figure  4  link  up  at  the  velocity  CL  -  4.51  km/sec  to  form  a  continuous 
curve  which  is  almost  flat  across  the  entire  frequency  range  of  0.02  to 
2.50  Hz.  The  separate  modes  thus  combine  to  form  a  single  non-disper- 
sive  phase  propagating  at  the  shear  wave  velocity  of  the  mantle  LVZ,  and 
an  examination  of  the  behavior  with  depth  of  the  particle  displacement 
eigenfunctions  reveals  that  this  phase  is  in  fact  the  mantle  LVZ  channel 
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TABLE  T 

Canadian  Shield  Earth  ModeL 
(from  Rrune  and  Dorman,  1963). 


Layer 

i  ickness  (km) 

Depth  at 

Bo  t  tom  ( km ) 

P-Wave 

VeLocity  (km /sec) 

S-Wa ve 

Velocity  km/sec ) 

Density 
(  pm/cm'  ) 

6.0 

6  0 

5.64 

3.47 

2  .  70 

10 . 5 

16.5 

6.15 

3.64 

2 . 80 

18.7 

35.2 

6 , 60 

3.85 

2 . 85 

80.0 

115.2 

8.  10 

4.72 

3  .  30 

100.0 

215.2 

8.20 

4 . 54 

3.44 

100.0 

315.2 

8.30 

4.51 

3.53 

80.0 

395.2 

8.70 

4 . 76 

3.60 

oo 

oo 

9  .  30 

5.12 

3.76 

SHERR  VELOCITY  C KM/SEC  I 
CRNRDIRN  SHIELD.  NO  SEDIMENT 


Figure  1,  Shear  wave  velocity  versus  depth  for  the  Canadian 
shield  model  with  no  sediments. 


CANADIAN  SHIELD.  NO  SEDIMENT 


OENSITY  C  GM/CM»*3  I 


Figure  2. 


Density  versus  depth  for  the  Canadian  shield 
with  no  sediments. 
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FREQUENCY  C  HZ  3 


CANADIAN  SHIELD.  NO  SEDIMENT 


Figure  4.  Phase  velocity  versus  frequency  for  twenty-six  Love 
wave  modes  of  the  Canadian  shield  model  with  no 
sediments . 


WaVG  ^n*  generation  of  by  the  superposition  of  higer  Love-wave 

modes  is  examined  in  detail  by  Stephens  and  Tsacks  (1977)  and  by 
Mantovani  et  al.  (1977).  Because  S  is  a  guided  wave  in  the  LVZ,  it  Is 
inefficient  at  exciting  motion  at  the  earth's  surface,  a  fact  which  is 
demonstrated  by  the  very  small  values  which  are  computed  for  the  propa¬ 
gation  medium  transfer  function  for  frequencies  at  which  the  Love-wave 
phase  velocities  are  greater  than  4.51  km/sec.  We  are  therefore 
Justified  in  limiting  our  synthetic  seismograms  to  the  velocity  range 

shown  in  Figure  4,  C  <  4.51  km/sec,  effectively  defining  "L  "  for  our 

..  g 

purposes  as  being  "that  portion  of  the  transverse  component  of  motion 

which  occurs  after  the  arrival."  The  cut-off  in  frequency  of  2.5  Hz 

was  chosen  for  computational  efficiency;  our  synthetic  seismograms  will 

have  this  cut-off  as  their  Nyquist  frequency,  and  hence  they  will  be 

generated  at  a  sampling  rate  of  5  samples  per  second.  This  data  density 

is  high  enough  for  a  resonable  comparison  with  L  as  observed  on  most 

ft 

short-period  seismographs.  Twenty-six  modes  were  treated  because  these 
are  all  the  Love-wave  modes  there  are  (for  this  particular  model)  for 
which  the  roots  of  the  period  equation  lie  within  the  two  ranges  < 

4.51  km/sec  and  f  <  2.50  Hz. 

It  is  instructive  to  compare  the  graph  of  the  phase 
velocity-versus-f requency  dispersion  relation  shown  in  Figure  4  with  the 
velocity  profile  listed  in  Table  I.  The  figure  shows  that  at  high 
frequencies  the  phase  velocity  of  the  fundamental  mode  asymptotically 
approaches  the  shear-wave  velocity  at  the  earth's  surface.  Since  the 
graph  is  almost  horizontal  at  this  velocity,  the  fundamental  mode 
becomes  almost  nondispersive ,  so  we  expect  that  the  energy  contained  in 
this  mode  will  arrive  in  almost  a  single  pulse  at  high  frequencies.  In 
several  figures  to  be  shown  later,  we  shall  demonstrate  that  the 
fundamental  mode  does  in  fact  have  a  prominent  Airy  phase  with  a  group 
velocity  of  3.47  km/sec  which  causes  the  synthetic  seismograms  for  this 
mode  to  be  impulse-like. 

Figure  4  shows  further  that  the  higher  modes  approach,  but  do  not 
reach,  osculatory  behavior  at  phase  velocities  equal  to  the  shear-wave 
velocities  in  the  second  and  third  layers  within  the  crust  (cf.  Table 
I).  This  nearly  osculating  behavior  reflects  a  tendency  of  the  sharp 


interfaces  between  the  crustal  layers  to  act  as  waveguides,  trapping  by 
I  successive  reflections  the  higher-mode  energy  within  the  individual 

layers  as  crustal  channel  waves.  Because  the  modes  do  not  actually 
osculate,  these  crustal  channel  waves  do  not  fully  develop;  we  shall 
3ee,  however,  that  the  medium  transfer  function  does  exhibit  minima  at 
the  frequencies  and  mode  numbers  corresponding  to  the  nearly  horizontal 
segments  of  the  graphs  in  Figure  4,  showing  that  these  points  of  near- 
osculatlon  represent  values  at  which  L  Is  inefficient  at  generating 
particle  motion  at  the  earth's  surface.  We  note  that  the  modes  would 
truly  osculate  and  crustal  channel  waves  would  be  formed  in  a  structure 
with  a  crustal  LVZ,  such  as  the  Basin  and  Range  province  of  the  south¬ 
western  United  States. 

Figure  5  shows  the  group  velocity  as  a  function  of  frequency  for 
the  twenty-six  Love-wave  modes.  Although  the  group-velocity  values  of 
course  represent  the  derivative  2"-df/dk  which  would  be  calculated 
numerically  from  the  values  of  ^(f)  «*  2^/1^  shown  in  Figure  4,  the 
group  velocities  were  in  fact  calculated  directly  from  Haskell's  (1953) 
matrices  using  Harkrider's  (1970)  computer  program.  The  curves  shown  in 
Figure  5  are  thus  an  analytic,  rather  than  a  numerical,  result  (although 
they  are  subject  to  the  same  numerical  iraprecisions  as  are  the  other 
quantities  which  are  computed  from  the  products  of  the  propagator 
matrices).  Inspection  of  the  figure  shows  that,  as  was  anticipated  on 
the  basis  of  Figure  4,  the  fundamental  mode  has  a  very  broad 
group-velocity  minimum  at  UL  **  3.47  km/sec.  The  higher  modes  have  a 
complex  pattern  of  stationary  points  at  which  Airy  phases  will  occur. 
Since  most  of  the  energy  will  propagate  at  the  velocities  corresponding 
to  the  Airy  phases,  Figure  5  shows  that  most  of  the  energy  in  the  L 
phase  travels  at  velocities  of  between  3.41  and  3.81  km/sec  in  the  hard 
rock  Canadian  shield  model.  This  range  exceeds  the  velocity  which  is 
commonly  taken  for  L  propagation  UT  -  3.54  km/sec.  We  shall  later 
adjust  our  earth  model  in  order  to  reconcile  this  discrepancy,  but  first 
we  shall  compute  the  synthetic  seismograms  and  thereupon  see  how  much 
energy  is  predicted  to  arrive  with  these  high  velocities. 
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Figure  5,  Group  velocity  versus  frequency  for  twenty-six  Love 
wave  modes  of  the  Canadian  shield  model  with  no 
s  ed imen  t s , 
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Some  idea  of  the  frequency-  and  modal-dependent  distribution  of 
energy  in  Is  given  by  the  medium  amplitude  response,  or  tranfer 
function,  which  is  shown  In  Figure  6.  The  transfer  function  is  a 
measure  of  the  efficiency  of  each  mode  for  producing  motion  at  the 
earth's  surface  as  a  function  of  frequency.  More  exactly,  treating  the 
layered  earth  structure  as  a  linear  network,  the  transfer  function  is 
the  frequency-dependent  ratio  of  an  output  (in  microns  of  particle 
motion)  at  the  surface,  to  input  (in  dynes  of  the  impulse)  which 
generated  that  motion.  We.  emphasize  that  this  is  not  equivalent  to  the 
surface  motion  which  would  be  generated  by  an  arbitrary  seismic  source, 
since  the  transfer  function  must  be  multiplied  in  the  frequency  domain 
by  an  appropriate  source  function.  For  a  double— couple  mechanism,  for 
example,  the  spatial  orientation  of  the  double  couple  determines  the 
relative  excitation  of  the  stress  and  particle  displacement  eigen¬ 
functions,  so  the  transfer  function  must  be  multiplied  by  a  quantity 
which  is  dependent  upon  frequency,  mode  number,  and  focal  depth.  What 
the  transfer  function  does  show  is  the  surface  motion's  dependence  upon 
frequency  and  mode  number  when  the  other  parameters  are  held  constant. 

In  particular,  it  shows  the  surface  response  to  an  event  which  is  itself 
at  the  surface,  since  in  this  special  case  the  stress  eigenfunctions 
vanish  and  the  (normalized)  particle  displacement  eigenfunctions  yield 
a  value  of  unity,  independent  of  frequency  and  mode  number.  This 
significance  of  the  transfer  function  is  reinforced  by  an  examination  of 
Figure  6,  which  in  light  of  the  previous  discussion  may  be  interpreted 
as  showing  that  for  a  surface  source  only  the  first  three  modes  make  a 
major  contribution  to  the  surface  motion  of  Lg  within  the  frequency  band 
of  interest.  This  is  hardly  surprising,  since  the  propagation  of  these 
modes  Is  confined  to  shallow  depths,  as  we  shall  subsequently 
demonstrate,  and  hence  they  are  more  easily  excited  by  a  surface  source 
than  are  the  higher  modes  which  propagate  at  greater  depths.  In  this 
discussion  we  are  using  the  "depth  of  propagation”  in  a  rather  nebulous 
sense,  which  we  shall  refine  later,  to  indicate  the  distribution  with 
depth  of  the  kinetic  energy  density.  Figure  6  also  shows  that,  all 
other  factors  being  held  constant  (which,  we  emphasize,  they  actually 

will  not  be),  most  of  the  low-frequency  energy  present  in  L  would  he 
.  ,  8 
contained  within  the  fundamental  and  first  higher  Love-wave  modes.  This 
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concentration  of  low-frequency  energy  will  be  especially  significant 
later,  when  we  alter  our  hard-rock  model  by  adding  a  sedimentary  layer 
through  which  these  shallow  modes  will  propagate.  Finally,  we  note  the 
pattern  in  Figure  6  of  prominent  minima  in  the  transfer  functions.  A 
comparison  with  Figure  4  reveals  that,  as  we  had  anticipated,  the  r ulls 
in  the  amplitude  response  occur  at  the  same  frequencies  and  mode  numbers 
(e.g.,  mode  7  at  1.56  Hz)  as  do  the  nearly  horizontal  segments  of  the 
phase  velocity-versus-f requency  graphs.  The  flat  portions  of  these 
graphs  correspond  to  channel  wave-like  behavior  of  those  modes,  which  at 
those  frequencies  propagate  mainly  in  waveguides  within  the  crust  and 
hence  are  inefficient  at  generating  motion  at  the  surface. 


In  addition  to  computing  synthetic  L  seismograms  using  the  earth 

8 

model  whose  parameters  are  illustrated  in  Figures  1-6,  we  will  wish  for 
the  sake  of  comparison  to  generate  L  synthetics  using  an  alternate 

O 

earth  model.  In  particular,  we  will  wish  to  examine  the  effect  upon  L 


propagation  of  a  sedimentary  layer  overlying  the  shield  structure.  This 

alternate  earth  model  of  a  tectonic  shield  with  a  sedimentary  cover  will 

be  applicable,  in  at  least  a  gross  sense,  to  many  regions  of  the  world, 

such  as  central  North  America,  and  we  anticipate  that  a  wide  range  of 

potential  characteristics  which  L  may  exhibit  will  be  bounded  by  the 

two  extremes  of  the  hard-rock  and  the  soft-rock  cases.  We  shall  thus 

consider  the  effect  of  replacing  the  topmost  1  km  of  the  hard-rock  earth 

model  described  in  Table  I  by  a  1  km-thick  sedimentary  layer  which  we 

shall  take  to  have  a  shear-wave  velocity  of  2.11  km/sec  and  a  density  of 

3 

2.39  gm/cm  .  These  values  were  deliberately  chosen  to  be  significantly 

smaller  than  those  in  the  hard-rock  case,  so  that  the  effect  upon  L  of 

g 

adding  the  sediments  would  be  prominent.  Figure  7  shows  that  adding  the 
sedimentary  layer  affects  the  phase  velocity  of  the  lowest-order  modes 
most  strongly.  This  behavior  of  the  dispersion  relation  is  predictable, 
since  it  is  the  lowest-order  modes  which  propagate  largely  in  the 
topmost  layer  and  whose  phase  velocity  asymptotically  approaches  at  high 
frequency  the  shear-wave  velocity  of  the  sediments.  The  higher-order 
modes  penetrate  deeper  into  the  crust,  so  they  are  less  perturbed  by  the 
presence  of  the  sedimentary  veneer. 
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Figure  7.  Phase  velocity  versus  frequency  for  twenty-eight 
Love-wave  modes  of  the  Canadian  shield  model  with 
1  km  of  sediments. 


Six 


The  sedimentary  layer  has  a  dramatic  effect  not  only  on  the  phase 
velocity  but  also  on  the  group  velocity  of  the  lowest-order  modes,  as  is 
shown  in  Figure  8.  This  graph  indicates  that  L  propagation  in  this 

s 

earth  model  will  be  characterized  by  two  strong  phases  arriving 

significantly  later  than  the  rest  of  the  energy  in  the  signal:  first  a 

low-frequency  (  ;  0.6  Hz)  Airy  phase  due  to  the  fundamental  mode,  and 
then  a  high-frequency  one  (  =  2.15  Hz)  due  to  the  first  higher  mode.  In 
the  discussion  which  follows,  we  shall  discuss  why  waveforms  with  this 
unusual  signature  will  noc  be  observed  in  practice. 

Figure  9  shows  that  not  only  do  the  two  Airy  phases  arrive  much 

later  than  the  rest  of  the  signal  but  they  also  are  much  greater  in 

amplitude  than  the  preceding  L  arrivals.  Specifically,  the  figure 
indicates  that  the  layered  medium's  amplitude  response  is  quite  large 
for  the  two  lowest  modes,  especially  at  frequencies  close  to  or  greater 
than  those  corresponding  to  the  group-velocity  minima  in  Figure  8.  We 
once  again  point  out  that  in  the  process  of  generating  the  synthetic 
seismograms,  the  values  shown  in  Figure  9  must  be  multiplied  by 
frequency-,  modal-,  and  depth-dependent  source  terms,  so  the  two 
prominent  lowest  order  modes  may  not  be  excited  for  a  given  particular 
case.  Figure  9  does  show,  however,  that  these  two  modes  will  be  more 
prominent  for  a  seismogram  generated  using  the  soft-rock  model  than  for 
one  with  the  same  source  mechanism  and  focal  depth  In  the  hard-rock 
model.  This  difference  should  especially  be  evident  for  shallow  events, 
which  excite  the  lowest  modes  preferentially,  and  in  particular  It 
should  be  prominent  for  events  with  hypocenters  within  the  sedimentary 
layer. 

In  Figures  10-12  the  effect  is  shown  of  increasing  the  thickness  of 
the  sedimentary  layer  from  1  km  to  2  km.  The  phenomena  which  were 
described  in  Figures  7-9  are  seen  to  be  enhanced,  but  the  changes  are 
less  dramatic  than  those  which  were  introduced  by  the  addition  of  the 
1  km  cf  sediments  to  the  hard-rock  model-  In  our  synthetics  we  shall 
use  the  2-km  thick  sediments  rather  than  the  1-km  thick  layer  in  order 
to  maximize  the  differences  between  the  hard-rock  and  the  soft- rock 
cases . 
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Figure  9.  Layered  medium  amplitude  response  for  twenty-eight 
Love-wave  modes  of  the  Canadian  shield  model  with 
1  km  of  sediments.  The  arrows  indicate  the  strong 
amplitude  response  of  the  Airy  phases  of  the  first  two 
modes,  which  correspond  to  the  broad  group-velocity 
minima  in  Figure  8. 
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Figure  10.  Phase  velocity  versus  frequency  for  twenty-nine 
Love-wave  modes  of  the  Canadian  shield  model 
with  2  km  of  sediments. 
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Figure  11.  Group  velocity  versus  frequency  for  twenty-nine 
Love-wave  modes  of  the  Canadian  shield  model 
with  2  km  of  sediments. 
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Figure  12.  Layered  medium  amplitude  response  for  twenty-nine 
Love-wave  modes  of  the  Canadian  shield  model  with 
2  km  of  sediments. 


Figure  13  Illustrates  the  effect  of  anelastic  attenuation,  which 
was  not  shown  for  the  hard-rock  model.  The  reason  for  ignoring 
attenuation  in  the  discussion  of  the  hard-rock  model  is  that  in  that 
particularly  simple  case,  for  which  the  shear-wave  quality  factor  Q  ^  was 
assigned  the  frequency-independent  value  1000  in  every  layer,  the 
trivial  result  is  obtained  that  •  1000  for  all  Love-wave  modes  and 
all  frequencies.  In  the  soft-rock  model,  however,  is  a  function  of 
depth,  and  is  a  complicated  function  of  frequency  and  of  mode  number. 
In  order  to  emphasize  the  effect  of  attenuation  in  the  soft-rock  case, 
we  have  set  Q  g  ■  20  in  the  2-km  thick  sedimentary  layer,  independent  of 
frequency.  This  anomalously  low  value  of  Qg  was  deliberately  chosen  in 
order  to  model  a  "worst  possible  case"  of  attenuation  by  a  sedimentary 
layer.  As  we  have  already  discussed,  the  use  of  the  hard-rock  model  in 
the  synthetic  seismograms  leads  to  the  appearance  of  prominent 
fundamental-mode  Airy  phases  which  are  not  observed  in  actual 
seismograms  of  L  .  One  way  to  diminish  the  contribution  of  the 

o 

fundamental  mode  to  the  synthetics  is  to  make  the  value  of  Q  much  lower 

in  the  topmost  layer  than  in  the  underlying  crust,  thereby 

preferentially  damping  out  the  shallowest  modes.  This  "worst  possible 

case"  model  will  serve  to  illustrate  the  effects  of  the  sedimentary 

layer  by  maximizing  them;  in  a  realistic  earth  model,  the  attenuation 

will  be  significantly  less  than  that  which  is  shown  in  Figure  13,  so  we 

expect  realistic  L  behavior  to  be  bounded  by  the  extremes  of  our 

hard-rock  and  soft-rock  limits.  As  typical  of  the  values  of  the 

anelasticity  quality  factor  which  one  might  expect  to  observe  in  the 

earth,  we  note  that  Bache  et  alj_  (1978)  inverted  Rayleigh-wave  data 

(which,  like  Love  waves,  are  sensitive  mainly  to  Q g  rather  than  ) 

measured  at  periods  of  T  >^3.0  sec  to  obtain  an  earth  model  for  the 

Nevada  Test  Site-to-Albuquerque  profile  consisting  of  a  2.5-kra  thick 

layer  for  which  Q  g-  20,  a  21.5-km  thick  layer  for  which  Qg  -  300,  and 

an  18-km  thick  layer  for  which  Qg  -  1000.  Since  we  are  concerned  with 

shorter-period  surface  waves  than  those  which  were  analyzed  by  Bache  et 

a1-  (1978),  we  expect  the  upper-layer  value  of  to  be  significantly 

higher  than  the  value  of  Q  -  20  which  they  found  and  which  is  used  in 

Figure  13.  We  shall  therefore  examine  alternative  methods  to  the 

introduction  of  an  anomalously  low  Q  in  order  to  diminish  the 

contribution  to  L  of  the  fundamental  mode, 

g 
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Figure  13.  Quality  factor  versus  frequency  for  twenty-nine 

Love-wave  modes  of  the  Canadian  shield  model  with 
2  km  of  sediments. 
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For  our  worst-case  (i.e.,  most  strongly  attenuating)  model,  we  set 

Qg"  20  Independent  of  frequency  In  the  surflclal  layer,  and  we  obtain 

the  frequency-  and  modal-dependent  values  for  Q  which  are  shown  In 

Li 

Figure  13.  These  curves  were  calculated  using  the  formula 

Q  l  "i  =  i  Pj  U^/C^)(8d^/3 Q  ^  for  N  layers  (1) 

(Knopoff,  1964).  The  behavior  of  this  function  Is  difficult  to  estimate 

Intuitively,  but  we  are  able  to  anticipate  certain  gross  features  of 

Figure  13  by  considering  the  behavior  with  depth  of  the  particle 

displacement  eigenfunctions  for  the  various  modes.  At  low  frequencies, 

the  fundamental  mode  sees"  (i.e.,  has  a  significant  particle 

displacement)  down  Into  the  upper  mantle,  and  so  It  is  not  greatly 

affected  by  the  attenuation  in  the  top  2  km  of  the  crust.  At  higher 

frequencies,  however,  there  is  little  fundamental— mode  particle 

displacement  at  depth,  and  so  the  mode  is  largely  confined  to  the 

attenuating  sedimentary  layer.  The  quality  factor  thus  effectively 

becomes  Q  g8ediment8,  as  is  shown  in  Figure  13.  A  similar  behavior  Is 

exhibited  by  the  first  higher  mode.  Figure  13  shows  that  Q  for  all 

L 

modes  of  number  3  or  higher  oscillates  with  increasing  frequency  before 
decaying  asymptotically  to  Q 6 sediment8*  The  reason  for  this  oscil¬ 
lating  behavior  is  that  the  eigenfunctions  for  these  modes,  which  will 
be  Illustrated  shortly,  oscillate  with  depth,  and  as  the  frequency 
increases  there  will  be  an  increasing  amount  of  a  vertical  "wavelength” 
(In  a  rather  loose  sense)  of  these  oscillations  contained  within  the 
topmost  2  km  of  the  crust.  If,  for  a  given  mode  number  and  a  given 
frequency,  the  particle  displacement  eigenfunction  Is  large  in  the  top  2 
km,  much  of  the  energy  of  that  mode  will  be  attenuated.  If  the 
frequency  is  increased  so  that  the  "wavelength"  of  the  oscillations  with 
depth  is  diminished,  the  eigenfunction  will  have  a  null  rather  than  a 
maximum  in  the  topmost  2  km,  so  most  of  the  energy  in  the  mode  at  that 
frequency  is  localized  below  the  attenuating  layer  and  Q  is  therefore 

ll 

'more  representative  of  the  lower  crust.  The  figure  of  merit  Q  thus 

Ij 

oscillates  in  frequency  as  alternating  maxima  and  nulls  in  the  eigen¬ 
functions  move  into  the  sedimentary  layer. 


Figure  13  shows  an  important  frequency  dependence  of  for  the 
lowest-order  modes  in  this  earth  model:  the  sedimentary  layer  causes 
the  crust  to  behave  like  a  low-pass  filter.  This  situation  contrasts 
with  that  of  the  higher-order  modes,  for  whic  is  not  a  monotonic 
function  of  frequency.  It  is  instructive  to  compare  the  quality  factor 
in  Figure  13  with  the  transfer  function  in  Figure  12.  The  graphs  of  the 
lowest-order  modes  exhibit  complimentary  behavior  in  the  two  figures, 
with  becoming  small  at  those  frequencies  where  A^  becomes  large. 

This  means  that,  as  was  indicated  earlier,  the  prominent  Airy  phases 
which  would  be  anticipated  on  the  basis  of  Figures  11  and  12  will  in 
fact  not  be  seen,  because  they  will  damp  out  rapidly  with  increasing 
source-to-receiver  separation  as  the  energy  which  1b  confined  to  the 
sedimentary  layer  diminishes  dues  to  anelasticity. 

The  soft-rock  values  of  which  are  shown  in  Figure  13  and  the 
constant  hard-rock  value  -  1000  will  be  applied  to  the  synthetic 
seismograms  by  computing  an  attenuation  coefficient 

YL  (f)  -  ^f/(QL  (f)  UL  (f)). 

This  coefficient  will  be  applied  to  the  amplitudes  by  multiplying  each 
point  in  the  frequency-domain  seismogram  by  exp  (-  Y  ^  r )  where  r  is  the 
source-to-receiver  separation  in  km.  Even  though  anelasticity  affects 
not  only  the  amplitude  but  also  the  velocity  of  seismic  waves,  we  shall 
not  consider  the  influence  of  Q  g  on  the  phase  and  ,?roup  velocities  of 
Lg,  since  this  is  a  small  effect  (Schwab  and  Knopff,  1971).  As  evidence 
of  this  smallness  (for  realistic  values  of  Qg),  we  note  that  at  the 

frequencies  under  consideration  body  waves  may  effectively  be  treated  as 

nondispersive.  Another  approximation  which  we  have  made  with  regard  to 
attenuation  is  that  Q  g is  independent  of  frequency.  Across  the 
frequency  band  under  consideration,  this  approximation  is  probably  a 
poor  one.  We  shall  nevertheless  disregard  the  behavior  of  with 

frequency,  since  we  are  more  interested  in  its  behavior  with  depth.  The 

approximation  Qg(f)  -  constant  Is  more  nearly  valid  than  is  the 
deliberately  exaggerated  approximation  Qg(h)  «  20  +  880H(h  -2),  where  h 
denotes  depth  and  H  is  the  Heaviside  step  function.  Our  soft-rock  model 
is  thus  meant  to  be  more  nearly  a  limiting  case  than  a  realistic  model, 
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and  we  shall  not  Introduce  higher  order  corrections  into  it.  Instead  we 
shall  present  alternative  methods  to  this  model  for  diminishing  the 
influence  of  the  fundamental  mode  at  shallow  depths. 

The  particle  displacement  and  stress  eigenfunctions,  to  which 

reference  was  made  in  the  preceding  discussion,  are  shown  in  Figures  14 

and  15  for  the  hard-rock  model  and  in  Figures  16  and  17  for  the 

soft-rock  model.  In  each  figure  the  eigenfunctions  are  shown  for  three 

specific  frequencies  and  for  several  specific  mode  numbers.  (Refer  to 

Figures  4  and  10  in  order  to  see  how  many  modes  exist  with  C.  <  4.51 

L 

km/sec  at  each  of  these  three  frequencies.)  We  wish  to  point  out 
several  features  of  these  plots  and  explain  their  significance  for  L 

g 

propagation.  First,  It  is  Important  to  note  that  the  transverse 
displacement  eigenfunctions  are  normalized  to  a  value  of  unity  at  the 
surface  and  that  the  different  modes  are  not  plotted  to  the  same  scale. 
The  stress  eigenfunctions  vanish  at  the  surface,  and  in  some  cases  this 
means  that  the  stress  changes  rapidly  from  zero  at  zero  depth  to  the 
first  value  which  is  shown  on  the  plots,  at  a  depth  of  0.1  km.  This 
rapid  change  will  not  concern  us,  since  we  shall  consider  no  focal 
depths  of  less  than  0.5  km,  and  it  is  the  value  of  the  eigenfunctions 
only  at  the  depth  of  the  hypocenter  which  enters  into  the  source 
function.  Next  we  note  that  the  eigenfunctions  are  in  certain  cases 
somewhat  ill  beha.jd  in  the  layer  immediately  above  that  layer  in  which 
they  cease  to  oscillate  and  instead  decay  exponentially  (i.e.,  the  depth 
at  which  B  >  C^).  Part  of  this  behavior  is  due  to  the  familiar  problem 
of  cumulative  numerical  imprecision  in  the  multiplication  of  a  sequence 
of  propagator  matrices,  and  part  is  due  to  the  true  behavior  of  the 
eigenfunctions  in  the  region  which  marks  the  transition  from  sinusoidal 
to  exponential  behavior  (the  "turning  point"  of  the  WKBJ  approximation). 
We  also  point  out  that  the  first  derivative  with  respect  to  depth  of  the 
tangential  stress  eigenfunctions  is  discontinuous  at  the  interfaces 
between  the  layers  defined  in  Table  I;  this  discontinuity  is  due  to  the 
discontinuities  of  the  modulus  of  rigidity  as  a  function  of  depth 
(cf.  Figure  3).  Finally,  we  point  out  that  Figures  14  and  16  illustrate 
a  phenomenon  which  we  have  mentioned  several  times  previously,  namely 
the  concentration  of  energy  at  shallow  depths  for  the  lowest-order  modes 
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at  high  frequencies.  These  figures  also  show  that,  as  the  frequency 
increases,  the  maxima  and  nulls  of  the  higher-order  modes  move  into  the 
near-surface  layers.  In  the  soft-rock  model,  when  the  particle 
displacement  eigenfunction  has  a  concentration  of  energy  in  the 
sedimentary  layer,  most  of  the  energy  in  that  mode  is  attenuated,  giving 
rise  to  the  minima  in  the  graphs  of  QL  in  Figure  13.  Similarly,  if  the 
energy  is  concentrated  in  the  underlying  shield,  Q  has  a  maximum.  This 
is  well  illustrated  hy  the  behavior  of  mode  5  at  0.48  and  1.48  Hz. 

Although  this  report  is  intended  to  investigate  the  effect  of 

crustal  structure  on  the  propagation  of  L  in  shields,  we  now  digress 

briefly  in  order  to  compare  the  results  of  our  computations  with  those 

which  would  result  from  using  a  very  different  earth  model,  namely  that 

of  the  Basin-Range  province  in  the  southwestern  United  States.  The 

shear-wave  velocity  profile  for  this  earth  model,  which  we  take  from 

Keller  et  al.  (1976),  is  shown  in  Figure  18.  An  important  difference 

between  this  model  and  our  shield  model  is  the  presence  of  a  crustal 

LVZ.  As  we  have  suggested  in  our  previous  discussions,  the  LVZ  causes 

the  phase  velocity  dispersion  curves  in  Figure  19  to  osculate  at  the 

shear-wave  velocity  of  that  layer.  The  high-velocity  "lid”  traps  energy 

in  the  LVZ,  with  the  osculating  modes  in  Figure  19  forming  a  continuous 

link  and  creating  a  nondispersive  wave  in  that  layer.  This  linkage  of 

modes  is  also  illustrated  in  Figure  20,  which  shows  that  the  transfer 

function  for  mode  3,  for  example,  decreases  by  more  than  six  orders  of 

magnitude  at  0.8  Hz,  where  the  transfer  function  for  mode  2  increases, 

and  then  it  increases  again  at  1.6  Hz,  where  the  transfer  function  for 

mode  4  plummets.  This  trade-off  causes  these  curves  to  link  up, 

yielding  a  relatively  flat  graph  for  A^  (f)  for  the  channel  wave. 

Although  we  cannot  make  any  quantitative  estimates  without  actually 

computing  the  synthetic  seisomgrams  for  the  Basin-Range  model,  we 

conclude  that  this  trapping  of  energy  in  the  crustal  LVZ  may  be  one  of 

the  principal  differences  between  L  propagation  in  shield  regions  and 

8 

in  the  southwestern  United  States.  The  presence  of  the  LVZ  correlates 

with  high  heat  flow  and  is  probably  indicative  of  partial  melting,  so  we 

anticipate  that  the  energy  in  this  layer  will  be  dissipated  by 

attenuation,  perhaps  resulting  in  a  stronger  L  attenuation  than  in 

8 
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Figure  18.  Shear  wave  velocity  versus  depth  for  the 
Basin-Range  model. 
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Figure  19.  Phase  velocity  versus  frequency  for  four 
Love-wave  inodes  of  the  Basin-Range  model 
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Figure  20.  Layered  medium  amplitude  response  for  four  Love— wave 
modes  of  the  Basin-Range  model.  The  fundamental  mode 
is  labelled  mode  1,  with  modes  2-4  as  the  first  three 
higher  order  modes. 


shield  structures.  Schwab  and  Knopoff  (1971)  appealed  to  similar 

reasoning  in  order  to  explain  the  failure  of  L  (at  longer  periods, 

8 

which  sense  the  upper  mantle  structure)  to  propagate  through  oceanic 
paths.  Although  we  have  considered  only  Love-wave  modes,  we  note  that 
the  high-velocity  lid  and  the  LVZ  will  also  trap  energy  which  propagates 
as  Rayleigh-wave  modes,  and  so  P-SV-type  head  waves  will  be  guided  along 
the  interfaces.  These  phenomena  may  contribute  to  the  amplitude  of 
crustal  phases  such  as  which  are  more  prominent  in  the  southwestern 
United  States  than  in  shield  regions;  there  may  in  fact  be  a  sort  of 
trade-off  whereby  the  crustal  structure  is  responsbile  for  partitioning 
energy  between  P  and  L  .  We  suggest  that  in  a  future  study  both 

o  8 

Love-wave  and  Rayleigh-wave  synthetic  seismograms  be  generated  for  the 

Basin-Range  model  in  order  to  investigate  the  generation  and  propagation 

of  regional  phases  in  a  tectonically  active  province.  In  particular, 

the  relationship  between  P  and  L  needs  to  be  understood  if  these  two 

8  8 

phases  are  to  be  used  as  seismic  discriminants. 

Returning  now  to  our  Canadian  shield  models,  we  shall  generate  syn¬ 
thetic  seismograms  using  the  methodology  which  we  have  described 
previously.  We  shall  use  these  synthetics  to  investigate: 

1)  the  frequency-  and  modal-dependence  of  L  for  various 
source  mechanisms  and  focal  depths;  ® 

2)  the  decay  of  L  with  distance  for  propagation  in 
the  hard-rock  §nd  soft-rock  models; 

3)  the  phase-velocity  structure  of  L  ;  and 

8  * 

4)  differences  between  synthetic  and  observed  L 

seismograms.  ® 


RESULTS 


Hard  Rock  Model 

We  examine  first  the  results  of  applying  the  previously  described 
process  of  seismogram  synthesis  to  the  "hard-rock"  Canadian  shield  model 
illustrated  in  Figures  1-6.  Figure  21  shows  the  contribution  of 
selected  modes  to  the  transverse  seismogram  which  would  be  recorded  on 
an  LRSM  short-period  instrument  at  a  distance  of  500  kra  from  a  0.5-km 
deep  strike-slip  earthquake.  The  figure  illustrates  each  of  the  first 
five  modes,  the  superposition  of  these  five  modes,  the  tenth  mode,  the 
superposition  of  the  first  ten  modes,  the  fifteenth  mode,  etc.,  and 
finally,  the  superposition  of  all  26  modes.  We  see  that,  as  we  had 
anticipated,  for  this  shallow  source  almost  all  of  the  energy  is 
transmitted  in  the  lowest— order  modes.  The  first  two  modes  show  very 
strong  pulses  where  energy  distributed  throughout  broad  ranges  of 
frequencies  arrives  almost  simultaneously.  These  spikes  in  the  L 
envelope  correspond  to  the  flat  parts  of  the  group  velocity  dispersion 
curves  (cf.  Figure  5),  3.46  km/sec  for  mode  1  and  3.42  km/sec  for  mode 
2.  Although  a  group-velocity  gate  was  imposed  on  the  synthetic 
seismogram  in  order  to  prevent  "wrap-around”,  some  numerical  error 
appears  to  be  responsible  for  the  late-arriving  (U  <  3.29  km/sec) 
energy,  since  no  such  arrivals  are  predicted  in  Figure  5.  The  principle 
source  of  the  computational  error  is  aliasing,  which  becomes  problematic 
at  frequencies  higher  than  about  f  £  2.2  Hz.  Because  the  LRSM 
instrument  response  attains  its  maximum  value  at  f  -  2.7  Hz,  there  is 
substantial  energy  in  the  wavetrain  at  frequencies  only  slightly  less 
than  the  2.5-Hz  Nyquist  frequency,  and  the  stroboscopic  effect 
introduced  by  the  sparse  sampling  of  this  energy  produces  "beats"  which 
can  be  seen  as  regularly  spaced  low-amplitude  peaks  in  the  otherwise 
featureless  portions  of  the  seismogram  both  before  and  after  the  passage 
of  the  large  wavepacket  containing  the  true  signal.  This  numerical 
problem  can  be  eliminated  by  suppressing  the  high-frequency  portion  of 
the  spectrum  by  means  of  an  anti-aliasing  filter  or  by  increasing  the 
sampling  rate  (in  which  case  more  modes  would  be  required  to  extend  the 
synthesis  to  the  new  Nyquist  frequency). 


The  behavior  of  the  higher  modes  can  be  seen  better  by  moving  the 
source  to  greater  depth,  as  is  done  in  Figure  22.  We  see  that  the 
influence  of  the  first  two  modes  is  much  smaller  for  this  10-km  deep 
source  and  that  the  U  “  3.46  km/sec  pulse  is  in  fact  all  that  can  be 
se  .1  of  the  fundamental  mode  when  all  the  modes  are  plotted  on  the  same 
scale.  The  second  mode  is  seen  to  increase  in  amplitude  as  the 
frequency  of  the  dispersed  wavetrain  increases  to  about  F  -  1  Hz,  and 
then  it  remains  approximately  constant,  in  accordance  with  the  behavior 
°f  (f)  for  this  mode,  which  is  shown  as  the  left-most  dashed  curve  in 
Figure  6.  The  dispersed  wavetrain  for  this  mode  terminates  abruptly 
when  all  the  high  frequencies  arrive  simultaneously  at  U  -  3.42  km/sec. 
The  fourth  and  all  higher  modes  exhibit  reverse  dispersion,  with  the 
lowest  frequencies  arriving  last;  again,  this  follows  from  Figure  5. 

Note  that  group  velocities  of  the  maximum  amplitudes  of  the  higher  modes 
are  less  than  those  of  the  lower-order  modes;  this  is  the  reverse  of  the 
situation  for  the  phase  velocities,  which  are  faster  for  the 
(deeper-penetrating)  higher-order  modes  (cf.  Figure  4).  Moving  the 
source  still  further  down  to  a  focal  depth  of  20  km  (Figure  23)  results 
in  the  complete  vanishing  of  the  fundamental  mode  and  in  the  significant 
propagation  of  energy  by  even  the  highest-order  modes.  Comparing  the 
traces  for  mode  5  in  Figures  22  and  23,  we  see  that  the  frequency 
content  of  this  mode  is  lower  for  the  deeper  source.  The  low-pass 
filtering  effect  is  brought  about  by  the  near  vanishing  of  the  stress 
and  displacement  eigenfunctions  for  mode  5  at  20  km  depth  for 
frequencies  greater  than  1  Hz,  and  the  existence  of  these  eigenfunctions 
at  frequencies  less  than  1  Hz,  as  is  shown  in  Figure  14a-c  and  15a-c. 

The  sum  of  all  26  modes  is  marked  by  sharp  spikes,  as  are,  to  some 
extent,  almost  all  of  the  synthetic  L  siesmograras  which  we  have 

s 

generated.  We  note  that  the  envelopes  of  observed  L  phases  are 
expected  to  be  smoother  than  those  of  the  synthetics  on  account  of 
scattering,  which  broadens  out  the  sharp  pulse-like  arrivals.  In 
Figures  21-23  we  have  considered  only  the  simple  case  of  strike-slip 
motion  on  a  vertical  fault.  The  shape  of  the  envelope  will  of  course 
change  as  the  focal  mechanism  is  varied  (while  the  depth  is  held 
constant),  as  is  shown  in  Figure  24.  Note  that  the  wavetrain  of  these 
seismograms  is  more  dispersed,  since  the  source— to— receiver  separation 
is  now  taken  to  be  1000  km. 
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Figure  24.  Synthetic  seismogram  for  different  focal 

mechanisms  at  a  depth  of  10  km  a  the  hard- 
rock  model. 


The  decrease  In  the  maximum  amplitude  of  the  L  envelope  with 

R 

increasing  focal  depth  is  illustrated  in  Figure  25  for  3  different 
combinations  of  focal  mechanism  and  receiver  distance,  all  for  the  same 
value  of  the  seismic  moment.  It  is  shown  that,  for  the  most  part,  the 
deeper  sources  have  a  smaller  Lg  amplitude,  and  that  the  sharpest 
decrease  in  amplitude  occurs  at  depths  between  5  and  10  km.  We  should 
note  that,  as  Figures  21-24  indicate,  the  visual  smoothing  of  the 
synthetic  wavetrains  to  form  envelopes  around  the  wave  packets  is  a 
rather  subjective  procedure,  so  the  maximum  amplitudes  plotted  in  Figure 
25  are  subject  to  some  uncertainty. 

We  now  wish  to  measure  the  rate  of  decay  of  the  maximum  amplitude 
with  increasing  source-to-receiver  separation.  Nuttli  (1973)  assumed 
that  the  maximum  is  an  Airy  phase,  in  which  case  its  amplitude  may  be 
expressed  as 


A  -  KA  (sin  A)  ^ ^  exp(-yA)  (2) 

where  K  is  related  to  the  seismic  moment  of  the  event  and  to  its 

efficiency  of  excitation,  A  is  the  angular  distance  from  the  source 

to  the  receiver,  and  y  Is  a  coefficient  which  we  shall  determine 

empirically.  In  Figures  26a-g  we  show  the  envelope  maxima  of  synthetic 

seismograms  generated  using  several  different  focal  mechanisms  and 

depths  in  the  hard-rock  model.  In  each  case  synthetics  were  generated 

at  distances  of  300,  600,  900,  1200,  1500,  and  1800  kra,  and  the  best- 

fitting  value  of  y  was  found  by  the  technique  of  least  squares  applied 

to  equation  (2).  The  amplitude  decay  which  is  computed  using  this  value 

of  and  Nuttli ' s  (1973)  formula  is  plotted  for  the  range  100  km  £  A  <_ 

2000  km  along  with  the  measured  amplitudes  of  the  synthetic  seismograms 

in  Figures  26a-g.  These  plots  show  that  in  every  case  some  value  of  Y 

can  be  found  so  that  Nuttli' s  (1973)  Airy-phase  formula  yields  a  close 

approximation  to  the  measured  values  of  the  envelope  maxima,  even  though 

these  maxima  were  not  necessarily  measured  at  the  same  group  velocity 

for  each  seismogram.  It  should  be  noted  that  the  strong  nonlinearity  of 

these  best-f itting  curves  (on  a  log-log  scale)  except  over  short 

distance  ranges  precludes  describing  the  amplitude  decay  in  terras  of  a 
-  a 

simple  A  falloff,  as  is  normally  done  for  longer~period  surface  waves 
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Figure  25.  Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  source  at  various  depths 
in  the  hard-rock  model, 

Circles:  rake  *  0°,  dip:  90°,  strike  =  0°  , 

A  *  500  km. 

Triangles:  rake  *  40°,  dip  =  70°, 

strike  -  30°,  A  =*  500  km. 

Squares:  rake  *  0°,  dip  -  90°,  strike  -  0°, 

A  *  1000  km. 
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Figure  26a.  depth  *  0.5  km, rake  =  0°,  dip  =  90°: 

Y  *=  0 . 1229  deg  1 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  hard-rock  model  as 
observed  at  various  distances. 
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Figure  26c. 
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Amplitudes  of  envelope  maxima 
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Figure  26d.  depth  -  20  kin,  rake  =  0°,  dip  =  90°: 

Y  -0.1785  deg"1 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  hard-rock  model  as 
observed  at  various  distances, 
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Figure  26e. 


depth  ■  30  km,  rake  =  0°,  dip  =  90°: 
y  =0.1778  deg"1 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  hard-rock  model  as 
observed  at  various  distances. 
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Figure  26E.  depth  =  10  km,  rake  =  90°,  dip  =  90°: 


y  =  0.1698  deg  1 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  hard-rock  model  as 
observed  at  various  distances. 
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(conventionally,  a  -  1.66  for  T  =  20  sec).  If  we  do  attempt  to 
approximate  the  observed  amplitude  decay  curves  by  straight  lines  on 
these  log-log  plots,  It  appears  that  perhaps  the  best  fits  are  given  by 
slopes  of  -1  In  the  distance  range  600  km  <.  A  <_  900  km,  -2  In  the 
distance  range  900  km  <  A  <1200  km,  and  -3  In  the  range  1200  km  1  A 
1800  km.  For  the  0.5  km-deep  source,  better  fits  are  given  by  slopes  of 
-1  for  300  km  i  A  f.600  km  and  by  -2  for  600  km£  A  iggg  km. 

A  potentially  valuable  technique  for  distinguishing  earthquakes 
from  explosions  Is  the  measurement  of  the  phase  velocity  of  the  crustal 
surface  wave  Lg.  We  anticipate  that  since  deeper  sources  preferentially 
excite  the  higher  modes,  they  will  generate  signals  which  propagate  at  a 
higher  phase  velocity  than  do  the  signals  from  shallow  sources,  which 
are  enriched  In  lower-order  modes.  Such  an  effect  has  been  observed  at 
the  Yellowknife  array  by  Barley  (1978).  As  a  crude  rule  of  thumb,  then, 
an  Lg  signal  with  a  high  phase  velocity  may  be  considered  as  evidence 
that  the  seismic  source  Is  located  at  intermediate  or  large  depths 
within  the  crust,  and  hence  that  It  Is  an  earthquake  rather  than  an 
explosion.  In  order  to  test  this  rule  of  thumb,  we  have  simulated  the 
measurement  of  Lg  by  a  closely-spaced  linear  array.  We  have  generated 
synthetic  seismograms  for  nine  receivers  located  along  a  line  at 
distances  of  500.0,  501.3,  502.2,  503.0,  503.5,  504.0,  504.8,  505.7,  and 
507.0  km  from  the  source  and  also  for  nine  other  receivers  located  at 

distances  of  1000.0,  1001.3 .  1007.0  km.  Either  one  of  these  two 

linear  arrays  can  be  used  to  compute  the  phase  velocity  of  a  dispersed 
wave train  by  means  of  frequency-wavenumber  (F-K)  analysis  (Barker  et 
— — — '  1980b)*  We  use  the  two  arrays  in  order  to  see  whether  the  phase 
velocity  Is  affected  by  either  dispersion  or  attenuation  as  the  signal 
travels  over  longer  paths.  The  results  of  the  phase-velocity  measure¬ 
ments  for  synthetic  seismograms  representing  several  different  depths 
and  focal  mechanisms  are  shown  In  Table  II.  We  have,  in  some  cases, 
computed  the  phase  velocity  for  two  separate  0.4-Hz  wide  frequency 
bands,  representing  both  a  particular  frequency  band  which  was  chosen  to 
remain  fixed  for  comparing  different  seismograms  and  another  band  which 
windows  the  frequency  determined  by  an  F-statistic  test  to  contain  the 
strongest  signal  for  each  individual  case.  It  Is  important  to  note  that 
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a  128-point  (i.e.,  25. 6-sec)  time-domain  window  was  chosen  which 
Included  the  maximum-amplitude  segment  of  the  signal.  Moving  the  window 
forward  or  backward  in  time  will  change  the  value  which  is  measured  for 
the  phase  velocity  for  each  frequency  band,  as  we  shall  demonstrate 
later.  Table  II  shows  that  our  rule  of  thumb  about  higher  phase 
velocities  for  deeper  events  has  some  validity,  but  there  is  a  large 
scatter  in  the  measurements  which  tends  to  blur  out  any  depth  dependence 
C.  •  We  have  already  noted  that  the  phase  velocity  is  dependent  on 
both  the  time  and  frequency  windows  which  are  selected  for  analysis;  we 
observe  from  Table  II  that  it  is  also  dependent  on  the  distance  from  the 
source  to  the  receiver  and  on  the  focal  mechanism  of  the  source.  In 
particular,  we  note  that  for  a  consistently  chosen  time—  and  frequency- 
domain  window  at  a  fixed  distance  of  500  km  from  a  10-km  deep  source, 
the  phase  velocity  varies  between  3.743  km/sec  for  a  strike-slip  fault 
dipping  at  60°  and  4.412  km/sec  for  a  strike-slip  fault  (or  for  one  with 
equal  strike-slip  and  thrust  components)  dipping  at  45°.  Since  this 
variability  at  a  single  depth  is  half  as  large  as  the  variability  over 
all  depths  shown  in  Table  II,  we  see  that  L  phase  velocity  perhaps  does 
serve  as  a  seismic  discriminant,  but  only  a  weak  one.  It  is  also 
possible,  however,  that  our  5  samples/sec  synthetic  seismograms  are 
Inadequate  for  showing  the  power  of  this  discriminant,  as  we  shall 
discuss  later. 

Soft  Rock  Model 

We  shall  now  repeat  the  preceding  analysis,  this  time  using 
synthetic  seismograms  generated  using  the  "soft-rock"  model  of  a  shield 
structure  overlain  by  strongly  attenuating  sediments.  Figure  27  shows 
that  the  low-Q  layer  effectively  absorbs  all  energy  propagating  in  the 
fundamental  mode,  even  for  a  shallow  source  which  tends  to  excite  that 
mode  preferentially  (cf.  Figure  21).  It  is  also  shown  that,  as  was 
anticipated  on  the  basis  of  Figure  13,  this  earth  model  acts  as  a 
low-pass  filter,  especially  for  the  lower-order  modes.  This  filter 
effect  is  also  evident  in  Figures  28  and  29,  which  show  the  wavetrains 
for  a  10-km  deep  earthquake  as  recorded  at  distances  of  500  and  1000  km. 
We  see  that  the  high  frequencies  become  less  and  less  prominent  the 
further  the  wavetrain  propagates.  Figure  30  shows  the  envelope  maxima 
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DEPTH  (km) 

Figure  30.  Amplitudes  of  envelope  maxima  of  synthetic 

seismograms  for  sources  at  various  depths  in  the 
soft-rock  mode  1  , 

Circles:  rake  =  0°,  dip  =90°,  strike  =  0°, 

A  =  500  km. 

Triangles:  rake  =  40°,  dip  =  70°,  strike  =  30°, 

A  =  5  0  0  km . 

Squares:  rake  =  0°,  dip  =  90°,  strike  =  0°, 

A  "  1000  km. 
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for  soft-rock  synthetic  seismograms  computed  using  the  same  focal 
mechanisms,  distances,  and  depths  as  were  used  for  the  hard-rock  model 
in  Figure  25.  The  amplitudes  are  no  longer  lrsrge  for  shallow  events, 
since  the  lower-order  modes  are  strongly  attenuated.  Tt  should  be  noted 
that  the  amplitudes  in  Figure  30  are  smaller  for  sources  at  a  depth  or  3 
km  than  at  depths  of  1.5  and  5  km.  The  signals  from  5— km  deep  sources 
are  stronger  than  are  those  from  3— km  deep  sources  because  the  deeper 
events  are  more  effective  at  generating  higher-mode  energy,  which  is 
less  attenuated  at  distance  than  is  the  energy  in  the  lowest-order 
modes.  The  signals  from  1.5-km  deep  sources  are  stronger  than  are  those 
from  3-km  sources  because  they  are  trapped  by  the  big  increase  in 
shear-wave  velocity  at  the  base  of  the  2-km  thick  sedimentary  layer. 

This  is  especially  true  for  the  case  denoted  by  triangles  in  Figure  30, 
since  for  events  with  this  focal  mechanism  the  exc  tation  of  Love  waves 
is  more  strongly  affected  by  the  low  value  of  the  modulus  of  rigidity  at 
the  hypocenter  than  are  the  pure  strike-slip  events,  which  are  denoted 
in  Figure  30  by  circles  and  squares  (Harkrider,  1970). 

Figure  31  shows  that  the  attenuation  is  in  fact  unrealistically 
high:  Y  1 8  about  three  times  larger  than  is  commonly  observed.  This 
too-large  value  of  Y  is  of  course  brought  about  by  the  low  value  of  Qp, 
independent  of  frequency,  which  was  used  for  the  2-km  thick  sedimentary 
layer  overlying  the  shield  structure.  Table  III  shows  that  the  phase 
velocity  cannot  be  used  effectively  as  a  depth  discriminant  for  signals 
in  this  model,  since  the  mesured  values  of  depend  not  only  on  depth 
but  also  on  A  .  As  the  waveform  travels  from  A  -  500  km  to  A  ■  1000  km, 
the  lower-order  modes  vanish  and  so  what  remains  are  the  higher-order 
modes  for  which  is  greater.  The  phase  velocity  thus  increases  as  the 
amplitude  decreases  with  increasing  source-to-receiver  separation. 

Constant  £  Model 

Aa  we  have  stated,  we  feel  that  the  anomalously  large  value  which 
we  found  for  Y  was  brought  about  by  our  using  too  low  a  value  for  in 
the  2-km  thick  sedimentary  layer.  In  order  to  isolate  the  effects 
introduced  by  the  elastic  parameters  B  and  P  in  the  sediments  as  opposed 
to  the  effects  of  Qg.we  have  generated  another  set  of  synthetic 
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depth  =  3.0  km1.rake  -  0°,  dip  *  90°: 

Y  =  0.2955  deg 

Amplitudes  of  envelope  maxima  of  synthet 
seismograms  for  the  soft-rock  model  as 
observed  at  various  distances: 
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Figure  31c 
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depth  =  5,0  km^  rake  =  0°,  dip  =  90°: 

Y  =  0 , 304 1  deg 

Amplitudes  of  envelope  maxima  of  syntheti 
seismograms  for  the  soft-rock  model  as 
observed  at  various  distances: 


Figure  31d.  depth  *  iO.O  km.  rake  =  0°,  dip  =  90°: 

Y  -  0,2601  deg"1 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  soft-rock  model  as 
observed  at  various  distances: 
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Figure  3 1  f .  depth  ■  30,0  kmj  rake  *=  0°,  dip  =  90°: 
y  =  0 , 2514  deg 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  soft-rock  model  as 
observed  at  various  distances: 
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Figure  31 g, 


depth  =  10.0  km.  rake  =  45°,  dip  =  60°: 
y  =  0.2760  deg 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  soft-rock  model  as 
observed  at  various  distances: 


seismograms  for  the  soft:  rock  model,  this  time  keeping  Q  constant  at 
1000,  even  for  the  sedimentary  layer.  The  dependence  of  amplitude  on 
depth  In  this  model  is  shown  in  Figure  32;  note  that  the  focal  mechanism 
represented  by  the  triangles  is  different  from  the  one  In  Figures  23  and 
30.  As  one  would  expect,  the  amplitudes  for  the  deep  (h  •  10  km) 
sources  are  roughly  the  same  in  Figures  23  and  32,  since  the  two  models 
appear  almost  identical  to  the  deeper-penetrating  higher  modes.  This  is 
in  contrast  to  the  situation  shown  in  Figure  30,  where  the  sedimentary 
layer  exerts  a  major  Influence  on  the  amplitudes  of  the  wavetrains  from 
even  the  deepest  sources.  We  conclude  that  this  influence  is  due  almost 
entirely  to  the  Q^tand  not  to  8  and  P  ,  in  the  sedimentary  layer.  For 
shallow  sources,  however,  we  see  that  the  presence  of  the  sedimentary 
layer  reduces  the  amplitudes  from  those  of  the  hard-rock  model  even  if 
the  value  of  Qpis  left  at  1000.  This  reduction  is  due  to  the  absence  of 
the  fundamental  mode  from  the  portion  of  the  wavetrain  In  which  T,  is 

g 

measured.  Figure  11  shows  that  the  modal  superposition  which 
constitutes  L  propagates  at  about  3.3  km/sec,  but  the  Airy  phases  of 

o 

the  four  lowest  order  modes  propagate  at  velocities  of  less  than  2.0 

km/sec.  The  fundamental  mode  will  thus  arrive  much  later  (by  107  sec  at 

"  500  km)  than  does  the  phase  which  is  taken  to  be  L  ,  so  it  will 

contribute  little  to  the  measured  L  amplitude,  regardless  of  Qp-  Figure 

11  shows  that  the  portion  of  the  fundamental  mode  which  Is  contained 

within  the  velocity  window  is  at  frequencies  too  low  to  be  seen  by  a 

short-period  seismograph,  so  L  in  the  soft-rock  model  will  not  contain 

8 

the  large  pulse-like  fundamental  mode  arrivals  which  characterized  the 
wavetrains  of  shallow  sources  in  the  hard-rock  model. 

We  note  that  in  practice  no  late  arrivals  are  actually  observed 
which  can  be  identified  as  uncontarainated  fundamental  mode  wavetrains. 

It  was,  of  course,  in  order  to  explain  the  absence  of  the  fundamental 
mode  that  the  strong  attenuation  was  introduced  into  the  model.  We  see 
in  Figure  32  that  simple  dispersion,  rather  than  attenuation,  is 

sufficient  to  remove  the  fundamental  mode  from  L  in  the  soft-rock 

8 

model;  we  must  still,  however,  find  a  way  to  account  for  the  absence  on 
real  seismograms  of  the  predicted  late-arriving  fundamental  mode 
wavetrain.  Although  introducing  low  Qpinto  the  sedimentary  layer  did 
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Figure  32,  Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  sources  at  various  depths 
in  the  high-Q  soft-rock  model. 

Circles:  rake  =  0°  ,  dip  =  90°,  strike  =■  0° 

A  =  500  km. 

Triangles:  rake  =  -30°,  dip  -  70°,  strike 

A  =  500  km. 

Squares:  rake  «  0°,  dip  *=  90°,  strike  -  0° 

A  =■  1000  km. 
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suppress  this  later  phase,  we  have  seen  that  it  led  to  unacceptably 

large  values  of  Y  for  I,  .  We  therefore  feel  that  a  more  likely 

8 

explanation  of  the  absence  of  the  fundamental  mode  is  scattering.  It  is 

at  least  a  viable  hypothesis  that,  on  account  of  topography  and 

geological  heterogeneity,  the  upper  layer  of  the  crust  is  more  effective 

at  scattering  the  surface-wave  energy  trapped  within  it  than  are  the 

lower  crustal  layers.  If  this  be  true,  then  scattering  will  more 

strongly  dissipate  the  shallow  fundamental  mode  than  the  deeper  modes 

which  contribute  to  L  in  the  soft-rock  model.  Scattering  will  also 

8 

prolong  the  L  coda,  thus  obscuring  any  sharp  fundamental  mode  arrival. 

8 

We  therefore  believe  that  a  more  accurate  model  of  L  propagation  could 

8 

be  formulated  by  using  a  higher  value  of  Qp.in  the  sedimentary  layer,  say 

perhaps  -200,  and  introducing  a  mechanism  for  scattering  at  shallow 

depths.  Figure  33  shows  that  a  more  nearly  realistic,  but  still 

somewhat  too  large,  value  of  Y  would  be  obtained  in  the  soft-rock  model 

if  Q^were  larger.  We  note  that  phase  velocity  is  still  a  poor  depth 

discriminant  in  this  model  even  if  1000  in  the  sedimentary  layer, 

since  L  (defined  as  the  phase  which  propagates  at  about  3 . 5  km/sec) 

8 

contains  little  energy  from  the  lower-order,  shallow  modes,  and  thus  the 
phase  velocity  is  high,  regardless  of  source  depth. 

A  comparison  of  the  synthetic  L  seismograms  with  real  waveforms 
has  been  presented  elsewhere  (Tier  et  al . ,  1981).  The  synthetic 
seismograms  which  were  used  in  that  report  differ  somewhat  from  those 
which  have  been  presented  herein,  since  the  upper  (i.e.,  sedimentary) 
layers  of  the  crust  were  modified  to  represent  the  structure  for 
individual  stations.  It  was  found  that  the  envelopes  of  the  synthetic 
seismograms  decay  faster  than  do  those  of  the  real  waveforms,  i.e.,  the 
synthetics  do  not  contain  long  enough  codas.  As  we  have  stated,  we 
believe  that  in  order  to  reproduce  on  the  synthetics  the  prolonged 
envelopes  which  are  characteristic  of  real  data,  it  will  be  necessary  to 
take  account  of  "he  late-arriving  energy  which  arrives  at  the  receiver 
by  scattering  ratuer  than  by  travelling  along  the  direct  path.  The 
observed  coda  shapes  exhibit  significant  variations  even  among  stations 
close  to  each  other,  thereby  showing  that  localized  site  effects  can 
dominate  the  behavior  of  L  .  F.vidence  for  the  effects  on  L  of  strictly 
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Figure  33d 


depth  =  10,0  km.  rake  =  0°,  dip  =  90°: 

Y  -  0.1555  deg" 

Amplitudes  of  envelope  maxima  of  synthetic 
seismograms  for  the  high-Q  soft-rock  modpl  as 
observed  at  various  distances: 
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Figure  3  3  f .  depth  =  30.0  km 
Y  =  0.1632  deg 
Ampl  i  tudes  of  e 
se i smograms  for 
observed  at  war 


local  structure  at  Yucca  Flats  and  Pnhute  Mesa  Is  presented  by  barker  et 
-a^  ‘  (1980a).  Our  synthetic  seismograms  are  thus  inadequate  in  that  we 
have  assumed  a  single  earth  model  which  is  presumed  to  describe  the 
structure  at  the  source,  throughout  the  entire  propagation  path,  and  at 
the  receiver.  Obviously,  our  approximation  of  uniform  plane-parallel 
layers  is  inadequate  over  long  distances,  so  our  synthetics  will 
reproduce  only  the  gross  behavior  of  L  .  An  important  phenomenon  in  L 

8  g 

propagation  is  almost  certainly  the  conversion  of  energy  between  modes, 
which  is  caused  by  scattering  and  by  the  departure  of  the  structure 
along  the  propagation  path  from  the  assumed  uniform  model.  We  have  made 
no  attempt  to  account  for  even  the  conversion  among  Love-wave  modes, 
much  less  the  conversion  between  Love  and  Rayleigh  modes  which  takes 
place  in  a  three-dimensional  scattering  problem.  We  note  that  taking 
into  account  the  Love-wave  energy  which  has  propagated  over  a  large 
fraction  of  the  path  as  (slower)  Rayleigh-wave  modes  would  extend  the 
codas,  since  most  of  the  scattering  along  the  path  occurs  at  shallow 
depths  aand  hence  tends  to  excite  the  fundamental  (i.e.,  slowest)  mode 
preferentially.  The  scattered  fundamental  mode  will  of  course  arrive 
at  the  receiver  after  the  higher-order  modes  but  before  the  direct-path 
fundamental  mode  which  is  prominent  on  our  synthetics  as  a  separate 
arrival  in  the  soft-rock  model.  Interference  between  these  multiple 
arrivals  will  tend  to  obscure  the  appearance  of  the  fundamental  mode  as 
a  separate  phase. 

Deficiency  of  the  Model 

We  have  mentioned  that  it  is  difficult  to  use  L  phase  velocity  as 

a  depth  discriminant  since  the  value  which  is  measured  for  C  will  vary 

L 

temporally  as  the  Airy  phases  of  different  modes  arrive  at  different 
times  throughout  the  coda.  In  Table  IV  we  demonstrate  this  difficulty 
by  using  synthetics  which  were  generated  using  a  crustal  model 
apppropriate  for  WN-SD,  a  station  located  on  thick  sediments  atop  a 
shield-like  structure.  It  is  shown  that  the  phase  velocity  is  strongly 
influenced  by  the  particular  25. 6-sec  time  window  which  is  chosen  for 
making  the  measurement.  The  F-statistic  is  also  shown  to  vary 
temporally  throughout  the  coda,  presumably  depending  upon  whether  most 
of  the  energy  within  the  given  time  window  propage  _s  as  pulse-like  Airy 
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TABLK  IV 


Temporal  Variation  of  the  Men. sored  L  Phase  Velocity 

R 


0.8-1  .2 

11  z 

1.0-1 

.  4  Hz 

1  .  2-1  .  6 

Hz 

Wind  ow 

Start 

Phase 

Phase 

Phase 

epth 

Time  F- 

Statistic 

Velocity 

F-Statistic 

Velocity 

F-Statistic 

Ve l oc i ty 

( km) 

(sec)* 

(km/ sec) 

( km/ sec ) 

(km/sec) 

3.0 

21 

54.2 

4.06 

35.7 

4.04 

55  .  3 

3.91 

25 

57  .  9 

4.23 

31.7 

4.04 

50.0 

3.91 

29 

51.9 

4.23 

37.7 

4.  18 

39  .  1 

4.02 

33 

52  .  3 

4.23 

43 . 2 

4.04 

36.8 

4.02 

37 

19.9 

3.76 

53.8 

3  .  78 

60.2 

4 . 26 

41 

63.3 

4.84 

89  .  2 

4.33 

122.5 

4 . 69 

45 

61  .8 

4.84 

110.3 

4,18 

132.7 

4.69 

49 

23.0 

3.91 

138.6 

4.  18 

59.0 

4 . 69 

1.0 

21 

86.4 

3.91 

54.3 

4.04 

65.0 

4.40 

25 

104 . 5 

3.91 

52  .  3 

4.04 

66.5 

4.40 

29 

106.7 

3.91 

53.5 

4.04 

75.0 

4.40 

33 

90.5 

3.91 

63  .  7 

4.18 

77.5 

4.40 

37 

65.2 

3.91 

78.7 

4.3  3 

3  9.7 

4.40 

41 

97  .  2 

4.06 

77.1 

4.49 

77.5 

4 . 54 

45 

103.5 

4.06 

103.4 

4.49 

47.8 

4 . 54 

49 

76.3 

4.42 

52 . 8 

4.49 

62  .  3 

4 . 54 

1.0 

21 

33.9 

4.06 

49.3 

4.33 

23.8 

4.14 

25 

36.4 

4.06 

63  .  1 

4.33 

27 . 6 

4.40 

29 

37  .  2 

4.06 

63.9 

4.33 

27 . 9 

4.40 

33 

23  .  7 

4.06 

31  .  5 

4.33 

24 . 2 

4 . 54 

37 

43.4 

4.06 

67 . 7 

4.33 

46 . 3 

4 . 54 

41 

43.8 

4.06 

63.4 

4.33 

47.8 

4 . 54 

45 

36.4 

4.06 

62.5 

4.33 

64.9 

4 . 54 

49 

92.6 

4.23 

95.4 

4.66 

40.9 

4.68 

*  relative  to  an  arbitrary  reference  point  approximately  21  sec 
before  the  first  motion. 


phases  or  as  dispersed  wavetralns,  which  are  less  coherent  across  the 

7.0-km  long  array  (since  our  frequency  bands  have  finite  bandwidth,  and 

hence  dispersion  occurs  within  them).  We  note  that  our  assessment  of 

the  value  of  this  depth  discriminant  well  might  change  if  our 

seismograms  were  sampled  at  the  more  conventional  rate  of  20  samples  per 

second,  in  which  case  the  time  windows  would  he  only  6. A  sec  in  duration 

and  would  therefore  more  nearly  contain  only  single  Airy  phases.  These 

shorter  time  windows  would  thus  measure  better  the  modal  structure  of 

L  • 

R 

One  would  suspect  that  a  more  realistic  value  of  the  attenuation 
parameter  y  would  be  obtained  if  we  allow  the  shear-wave  quality  factor 
Q  3  to  vary  with  frequency.  In  a  study  involving  synthesis  of  the 
Rayleigh-wave  modes  of  L  ,  Bache  et  ah  (1980)  found  abnormally  large 

O 

values  of  y  for  a  model  representing  the  eastern  United  States,  a  result 
which  is  compatible  with  the  overestimation  of  y  which  we  found  for  the 
Love-wave  modes  in  the  low-Q  soft-rock  model.  Bache  et  al.  (19oG) 
speculated  that  introducing  a  frequency-dependent  Q^into  their  earth 
model  would  result  in  more  nearly  realistic  synthetic  seismograms,  but 
in  a  subsequent  study  (Bache  et  ah,  1981)  they  found  that  making  this 
change  resulted  in  synthetics  which  still  suffered  from  the  same 
inadequacies  as  before.  They  concluded  that  scattering,  rather  than 
frequency-dependent  attenuation,  is  the  best  explanation  for  the 
observed  rate  of  decay  of  amplitude  with  distance.  As  we  have  stated, 
we  believe  the  neglect  of  scattering  to  be  a  principal  deficiency  of  the 
seismograms  synthesized  by  modal  superposition,  and  until  some  account 
is  made  for  this  deficiency,  we  are  reluctant  to  introduce  further 
refinements  into  the  process  such  as  increasing  the  sampling  rate  and 
the  number  of  modes,  adding  more  layers  to  the  earth  model,  or  making 
assumptions  about  the  variation  of  Q  ^  with  depth  and  with  frequency. 


CONCLUSIONS 


We  have  attempted  to  analyze  the  propagation  of  the  seismic  phase 
kg  or<^er  to  increase  our  understanding  of  how  Its  propagation  across 
a  tectonic  shield  affects  its  observed  characteristics  at  the  receiver. 
It  Is  necessary  to  understand  the  propagation  path  effects  in  order  to 
determine  what  information  can  be  determined  about  the  seismic  source 
from  observations  of  L  .  Our  synthetic  seismograms  have  involved 

o 

exclusively  the  Love-wave  modes  (l.e.,  the  SH  component)  of  L  ,  and  we 

8 

have  thus  been  unable  to  examine  seismograms  corresponding  to 
explosions,  since,  in  theory,  purely  radial  source  functions  produce  no 
transverse  motion.  Restricting  our  study  to  earthquakes,  we  have 
reached  certain  conclusions  about  the  propagation  of  L  across  shield 
structures : 

o  There  is  a  tendency  for  the  phase  velocity  to  increase  with 
depth,  but  the  potential  use  of  this  tendency  as  a  seismic 
discriminanc  to  identify  deep  events  (i.e.,  earthquakes)  is 
obscured  by  large  variations  in  the  phase  velocity  which  are 
due  to  focal  mechanisms.  The  measured  value  of  the  phase 
velocity  also  depends  on  the  exact  time  window  which  is 
processed  and,  in  a  structure  overlain  by  strongly  attenuating 
sediments,  on  the  distance  from  the  source  to  the  receiver. 

o  A  good  fit  to  the  observed  araplitude-versus-distance 

relationship  is  given  by  Nuttli's  (1973)  Airy-phase  formula, 
which  matches  the  amplitudes  of  the  synthetic  seismograms  over 
a  greater  distance  range  than  does  any  simple  A-0  falloff. 

In  order  to  use  Nuttli’s  (1973)  formula  to  calculate  the 
seismic  magnitude  tram  the  observed  amplitude,  however,  it  is 
necessary  to  evaluate  the  anelastic  attenuation  Y  ,  and  our 
attempts  to  estimate  this  parameter  in  an  earth  model  with  a 
strongly  attenuating  upper  crust  resulted  in  its 
over estimation .  Although  we  might  try  varying  the  assumed 
depth  dependence  of  Qg,  we  are  skeptical  of  being  able  to  do 
so  with  sufficient  accuracy  to  estimate  Y  adequately  for  any 
particular  propagation  path  for  which  this  parameter  has  not 
already  been  determined  by  direct  observation. 

o  Synthesis  of  L  by  modal  superposition  is  inadequate,  in  that 
discrepancies  Between  observed  and  predicted  envelope  shapes 
are  introduced  by  departures  of  the  actual  earth  structure,  on 
both  a  regional  and  a  purely  local  scale,  from  the  assumed 
laterally  homogeneous  model  consisting  of  plane-parallel 
layers.  Our  synthetics  predict  strong  excitation  of  the 
fundamental  mode,  which  will  appear  as  a  strong  pulse  for 
shallow  sources  in  a  shield  model  without  sediments  or  as  a 
later  arrival  in  a  model  with  a  sedimentary  veneer.  In  fact, 


this  strong  fundamental  mode  Is  not  observed,  but  Introducing 
a  thick,  low-Qp  layer  Into  our  crustal  model  In  order  to  damp 
It  out  leads  to  too  large  a  value  of  Y  .  We  conclude  that  a 
more  likely  explanation  for  our  failure  to  observe  the 
fundamental  mode  is  that  It  Is  scattered  more  strongly  by 
topography  and  by  shallow  lnhomogeneltles  than  are  the 
higher-order  modes.  The  effects  of  the  propagation  path 
structure  on  the  observed  characteristics  of  L  cannot  be 
fully  understood  until  the  scattering  among  mo§es  and  between 
Rayleigh-  and  Love-wave  modes  is  taken  Into  account. 

Seismogram  synthesis  which  Includes  scattering  would  attempt 
to  reproduce  observed  seismograms  in  a  stochastic,  rather  than 
a  deterministic,  sense. 
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Attn:  Mr.  J.  Murphy 
11800  Sunrise  Valley  Drive 
Suite  1212 
Reston,  VA  22091 

Teledyne  Geotech 

ATTN:  Dr.  Z.  Der  &  Mr.  W.  Rivers 
314  Montgomery  Street 
Alexandria,  VA  22314 

Woodward-Clyde  Consulants 
ATTN:  Dr.  Larry  Burdick 
556  El  Dorado  St 
Pasadena,  CA  91105 

Weidlinger  Associates 
ATTN:  Dr.  J.  Isenberg 
3000  Sand  Hill  Road 
Menlo  Park,  CA  94025 


NON-U. S. 

National  Defense  Research  Institute 
ATTN:  Dr.  Ola  Dahlman 
Stockholm  80,  Sweden 

Blacknest  Seismological  Center 

ATTN:  Mr.  Peter  Marshall 

Atomic  Weapons  Research  Establishment 

UK  Ministry  of  Defense 

Brimpton,  Reading  RG7-4RS 

United  Kingdom 

NTNF  N0RSAR 

ATTN:  Dr.  Frode  Ringdal 
P.0.  Box  51 
N-2007  Kjeller 
Norway 


To  be  determined  by  the  project  office 


